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Analogue models for gravity intend to provide a framework where matter and gravity, as well as 
their intertwined dynamics, emerge from degrees of freedom that have a priori nothing to do with 
what we call gravity or matter. Bose Einstein condensates (BEC) are a natural example of analogue 
model since one can identify matter propagating on a (pseudo-Riemannian) metric with collective 
excitations above the condensate of atoms. However, until now, a description of the "analogue 
gravitational dynamics" for such model was missing. We show here that in a BEC system with 
massive quasi-particles, the gravitational dynamics can be encoded in a modified (semi-classical) 
Poisson equation. In particular, gravity is of extreme short range (characterized by the healing 
length) and the cosmological constant appears from the non-condensed fraction of atoms in the 
quasi-particle vacuum. While some of these features make the analogue gravitational dynamics of 
our BEC system quite different from standard Newtonian gravity, we nonetheless show that it can 
be used to draw some interesting lessons about "emergent gravity" scenarios. 
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Analogue models for gravity have provided a powerful tool for testing (at least in principle) kinematical features 

fj ■ of classical and quantum field theories in curved spacetimes [l[. The typical setting is the one of sound waves 

^h| propagating in a perfect fluid [2,l3|. Under certain conditions, their equation can be put in the form of a Klein-Gordon 

I i equation for a massless particle in curved spacetime, whose geometry is specified by the acoustic metric. Among the 

(^ various condensed matter systems so far considered, Bose-Einstein condensate (BEC) [J, Q had in recent years a 

prominent role for their simplicity as well as for the high degree of sophistication achieved by current experiments. In 

a BEC system one can consider explicitly the quantum field theory of the quasi-particles (or phonons), the massless 

excitations over the condensate state, propagating over the condensate as the analogue of a quantum field theory of 

^^ ' a scalar field propagating over a curved effective spacetime described by the acoustic metric. It provides therefore a 

^-H , natural framework to explore different aspects of quantum field theory in various interesting curved backgrounds (for 

^\ . example quantum aspects of black hole physics 0, 01 or the analogue of the creation of cosmological perturbations 

^^ ' ISIj ISl)- Unfortunately, up to now, the analogy with gravity is only partial: there is no analogy with some sort of 

r~^ , (semiclassical) Einstein equations, since it has not been possible to put the fluid equations, which are those describing 

f^ ' the dynamics of the acoustic metric, in a geometrical form which could eventually lead to a complete dynamical 

00 , analogy with general relativity [IQ] . Our aim is to fill this gap and to gather from the description of the gravitational 

^D ■ dynamics general lessons about possible features of "emergent gravity" scenarios. 

^ ' In BEC, the effective emerging metric depends on the properties of the condensate wave-function. One can expect 

" '^ . therefore the gravitational degrees of freedom to be encoded in the variables describing the condensate wave-function 

rS I , which is solution of the well known Bogoliubov-de Gennes (BdG) equation. The dynamics of gravitational degrees 

2 . of freedom should then be inferred from this equation, which is essentially non-relativistic. The gravitodynamics of the 

- - ' BEC should therefore be the analogue of some sort of Newtonian gravity, and we shall reinterpret the BdG equation 

as a modified Poisson equation. 

The "emerging matter", the quasi-particles, in the standard BEC, are phonons, i.e. massless excitations. Since we 
expect the quasi-particles to be the matter source in the Poisson equation, we run a priori into a problem: massless 
particles are not treatable in the framework of Newtonian mechanics. To avoid this issue, we shall then introduce 
a new term in the BEC Hamiltonian which will softly break the usual U{1) symmetry and therefore will allow the 
quasi-particles to acquire mass. 

The plan of the paper is as follows. In the first section, we discuss the dynamics of the condensate and of the quasi- 
particles in the particular case of a BEC with a soft L''(l) symmetry breaking term. We show how this modification 
affects the dynamics of the condensate, leading to a generalized BdG equation, and how it provides a mass to the 
quasi-particles. In the second section, we focus on identifying the analogue of the gravitational potential and we show 
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that it has to be identified with inhomogeneities in the condensate density. We then show that the BdG equation 
can be cast in a form that closely mimic a Poisson equation. By looking at different sources we shall then identify 
the analogue of the Newton constant (and its possible dependence on the momentum of the quasi-particles) as well 
as the cosmological constant. The final outcome of our analysis is that the gravitodynamics of the BEG with massive 
quasi-particles is encoded in a modified Poisson equation characterized by two main features. First gravity is of very 
short range, i.e. a localized particle generates a Yukawa type potential. Secondly the cosmological constant naturally 
appears from the non-condensed fraction of atoms in the quasi-particle vacuum. 

I. A MODIFIED BEG 

Let us start recalling the standard setup for a BEG of many atoms in a box of volume V, in the dilute gas 
approximation. In this limit it is possible to describe the atoms via a second-quantized field operator 



k 

whose evolution is encoded in the Hamiltonian Hq [il| 

Ho^J *t(x) (^- Av^ - ^ + ||^|2^ ^(x)d3x, (2) 

which generates the standard non-linear equation^ 

zn^* = [iJo,*] = -|^V2^-M^ + K|*|2^. (3) 

The operator ^ has dimension L^^'^: the quantum average of its modulus square on a given state represents the 
number density of atoms. The mass m is the mass of the atoms. The energy fi is the chemical potential and has 
dimension ML'^T~^ . The constant k, of dimension ML^T'"^, represents the strength of the two-bodies interaction 
between atoms. We neglect higher order contributions to the second quantized Hamiltonian. 

While the condensation process can be easily understood as a macroscopic occupation number of an energy level, 
there are several approaches to describe it mathematically. The mean field approach is particularly convenient: we 
say that the system of N bosons has condensed whenever the field VP develops a non-zero vacuum expectation value 
(vev) 

(f]|*|f7)=V, (4) 

where ip is a, classical complex field describing the condensate and it is sometimes called the condensate wave- function. 
If this mean field is non-vanishing, we have that the two point correlation function 

G(x,y) = (*t(x)#(y)), (5) 

tends to a non-zero constant when x, y are infinitely separated, i.e. the system develops long range correlations |ll| . 

However, the mean field approximation is not the most rigorous method to treat a BEG system. In particular, the 
particle-number-conserving approach has been proven to give more accurate predictions for the physical properties of 
the condensate in settings in which the number of atoms, N, is fixed. This dimensionless quantity, then, is used to 
expand systematically the equation for the evolution of the operator 'J, in powers of N^'^. While this method was 
shown to provide a more accurate description of BEG systems, it has also shown that the mean field approximation 
gives already very good predictions for quasi-static configurations |12l | . 

In this paper, we consider an idealized case where the confining potential is almost constant, both in space and 
in time. Moreover, we neglect boundary effects due to the finite size of the trap, and we assume that all the 



^ The field operator 'I' is obeying the commutation relations: 

[*(x), *t(y)] = 53(x - y), [*(x), -^(y)] = [>&t(x), *t(y)] = 0. 



other experimental parameters are weakly time-dependent. In this case, then, the mean field approximation is well 
motivated. 

In order to derive the properties of the BEC, we split the field operator in the condensate part ip (the mean field) 
and an operator x representing the atoms which arc not in the condensate: 

* = # + X, (6) 

and use this form of the operator in the evolution equation ([3]) . 

For the reasons discussed in the Introduction and further analyzed in section |lll we need to slightly modify the 
standard Hamiltonian Hq by introducing a term which is (softly) breaking the U{1) symmetry in ([3]). 

Ho^H^Ho + H),, ffA--^ /"d3x(*(x)#(x) + *t(x)^t(x)). (7) 

The parameter A has the same dimension as ji. With this new Hamiltonian, the non-linear equation ^ becomes 

5 . , . S^ - - ~ ~ 

ih—^ = [H, *] = -7-— V^* - At* + k|*P* - A*^ (8) 

We shall show in section (|IB[) that H\ generates a mass for the quasi-particle^. Even though H\ both creates and 
destroys pairs of atoms, it is not difficult to check that H\ is not commuting with the number operator N, 

[Hx:N] =-A /"d^x('*(x)*(x)-^t(x)^t(x)^ (9) 

while unitarity is preserved. In fact, when applied on a state with a definite number of atoms n we have: 

\n)^\n^2) + \n + 2), (10) 

which means that an eigenstate of the number operator is taken into a superposition of states with different occupation 
numbers^. However, the expectation value of the number of operator on its eigenstates is still constant 

ih—{n\N\n) = {n\[N,H]\n) = (n\[N,Hx]\n) ex (?i|n - 2) - {n\n + 2) =0. (11) 

Finally, we want to discuss some plausible physical setup corresponding to the U{1) symmetry breaking in BECs. 
Given the relationship existing between the Hamiltonian ([7]) and the number operator, a natural situation would be 
represented by a condensate which is able to exchange particles with some sort of reservoir, in such a way to preserve, 
on average, their number. Several settings in this sense could be conceived, e.g. with coupling with suitably tuned 
lasers. Another interesting concrete example is represented by the case of magnon condensates, see [l^] and references 
therein. 

A. Condensate 

We consider the dynamics generated by ([5]), from which we want to extract the equation of motion for the condensate 
ijj. The evolution of the mean field -0 is easily determined in terms of the eigenstates \E) of the Hamiltonian H: 

i^i^ = i^ii{E\ME)) = {E\ih^ME) = -|-VV - /i^ - M* + A^?i> + ZkubV + KmsV*, (12) 

at at at 2m 

where ms = {E\x^\E),nE = (Elx'xlE) encode the effect of the non-condensate atoms. This is the generalization of 
the Bogoliubov-de Gennes (BdG) equation for the condensate wave-function to the case A ^ 0. 



^ If we were discussing a relativistic field theory this massive particle would be nothing else than a pseudo-Goldstone boson. 

^ As a consequence, the particle- number-approach cannot be used as it is defined in Il2ll . A more general treatment should be developed 
in order to include these interesting situations where the number of atoms is not a conserved operator. This extension could also provide 
a tost for the accuracy of the mean field approximation, and hence on the stability of the condensation, in the U(l) breaking case. 



If we have N particles '^ in the condensate, the number density of the non-condensate fraction is of order 1/A^ with 
respect to the number density of the condensate. In particular, the terms m, n are of order 1/iV. At zeroth order, we 
have the generalization of the Gross-Pitaevskii (GP) equation: 

ot 2m 

The time independent homogeneous solution to the GP equation is 

„, = |^|2^/i±^, (14) 

K 

where we have fixed the phase of the condensate to be zero^. 

We define the healing length ^ as the length scale at which the kinetic term is of the order of the self-interaction 
term: 

ft' .2 ft' 

This length represents the spatial scale needed for the condensate to pass from the value nc = at the boundary of 
the region where it is confined to the bulk value ric- In other words it represents the scale of the dynamical processes 
involving the deformation of the condensate wavefunction. 

B. Quasi-particles 

The equation of motion for the particles out of the condensate is obtained by subtracting the equation for the 
condensate ([T^ from the equation for ^ given in ([5]). We are interested in the propagating modes, so we neglect the 
self-interactions. We obtain: 

ift^X = "I^V^x + {2K\ij\' - m)x + M' - A)x^ (16) 

Let us consider the case of homogeneous condensate with density ric given above. In this situation we have: 

d h'^ 

^f'7r,x^-T^^^x + U^ + 2^)x + i^x^- (17) 

ot 2m 

If we decompose the field x in its plane wave components, we can rewrite this equation as 

d h^P 

*ft oT^k = -^ — ak + (M + 2A)ak + /iaLk- (18) 

The mixing between a and a^ due to the evolution in time becomes then apparent. We therefore pass to the quasi- 
particle operators 0(x) 



;^E^ke*-, (19) 



k 

which are related to the particle operators through the Bogoliubov transformation 

ak = a(k)Sk + /3(k)S_k, with a"^ (k) - (3'^ (k) ^ I . (20) 

a, (3 are only functions of A; = |k|, since the condensate is homogeneous and isotropic. The equation of evolution for 
the quasi-particles is then given by 

ih^k^£{k)W, (21) 



* Throughout the paper we shall intend by N the average number of particles which is conserved (cf. equation ((TTJ). 

^ In appendix |X] we show that this is not an arbitrary choice, but rather a consequence of the situation we want to describe. 



with the energy 

f(h\ - 



£{k) = (p^ + ^\{^Ji + x) + i^±^tee] . (22) 



The BogoHubov coefRcients are given by: 



_ A{k)+£{k) Q2,j^.^_ M^ .23) 

"^''^- 2E{k) ' '^ ^''^ - 2E{k) A{k) + E{ky ^^^> 



where we have introduced the quantity 



The high energy hmit of these coefficients is: 



21,2 



A{k) = —— + /i + 2A. (24) 

Zra 



hni a^{k) = 1, hm (3^{k) = 0, (25) 

k — *cxj k — >oo 

which means that at large wave-number (and hence large momentum), the quasi-particle operators coincide with 
the particle operators. This matches the behavior of the energy, which becomes just the energy of a non-relativistic 
particle of mass to, just like a free atom. The dispersion relation ([221) suggests the introduction of the following 
quantities: 

c^ = ^^±^, M^^A^^^^m\ (26) 

TO (/x + 2A)2 ^ ' 

Cs plays the role of the speed of sound, while Ai plays the role of a rest mass for the quasi-particle. Since M is 
proportional to A, we clearly see that it is the term H\ that generates the mass of the quasi-particle. When A ^ 0, 
that is when H —>■ Hq, the quasi-particle becomes massless, i.e. a phonon, and the speed of sound reduces to the usual 
one in BEC. 

Notice that, in order to have a non-negative mass square, we have to require A > 0. In standard BEC, we usually 
assume that the chemical potential fi is positive: indeed if it was not positive, there couldn't be any condensation. 
In our case, we can relax this requirement and obtain that /x > — A as a condition. In the following we consider 
/i > 0, in order to be able to consider the case in which the correction we are inserting is very small, without affecting 
dramatically the condensation. Indeed, it is easy to see that a condensation can take place in the system with a soft 
C/(l) breaking by checking the behavior of the two points correlation function G(x, y) — (^^^(x)^(y)), which is still 
showing the presence of long range correlations, since the mean field tp is non- vanishing (cf. equation P^ V 

A4 is proportional to to, the mass of the atoms. By introducing the ratio C = A//x, we introduce the function -F(C) 

M^ ^ F(C)to2 = 4i^i±^TO^ (27) 

Under our assumptions, we have that C ^ 0. It is then straightforward to check that on this domain F((^) is a 
monotonic (increasing) function and that 

F(0) = 0, lim F(C) = 1. 

We conclude therefore that the mass of the quasi-particles Ai is always bounded by the mass of the atoms, Ai £ [0, to). 
It is also interesting to notice that using the variable C, the speed of sound is: 

1 + 2C^^ (28) 

" l-^C TO ^ ^ 

For C, small, we then have cj. ~ Knc/m,, which is the standard result, while, for C — > cx), c^ -^ 2Knc/m. 



C. The various regimes for the MDR 

Before moving on to the gravitational dynamics, let us discuss briefly the content of the dispersion relation (|22p for 
the quasi-particles, rewritten using Cg and A4. 

m = {^+-y+M'-t)'\ (29) 

where we are using the obvious notation p — hk to simplify the shape of the expressions. Let us define the characteristic 
momenta pA, Pb and pc such that 

4 4 

Pa __ „2 2 Pb _ ka2a 



2 2 t-B _ A^2 4 2 2 _ A^^„4 

CsPAj 1 ^ — ■'^' Cj, CgP(j — JVi Cg 



4m^ " ' 4m^ 



so that they are explicitly 



pI = Am^cl pI = 2mMcl pi = M^cl- 
They are related through the relations 

pI = 2F{0p1^^f\0p\- 

If C <C 1, which will be the regime we shall consider, we have also that 

Pc<PB<PA- (30) 

Taking into account (|30p. the characteristic momenta define different regimes: 

• If p 3> PA, the term p^ dominates, the dispersion relation (j29|) is well approximated by £^ ^ p^ /2ra, we are in 
the transphononic regime. 

• If on the contrary we have pc <C p ^ pa, we can safely neglect the term of order p*, we are then in the relativistic 
regime since the dispersion relation (|29|) is well approximated by £ ^ (p^c^ + M^c^)^ ■ The quasi-particle is 
then relativistic, when the speed of sound Cg is playing the role of the speed of light. 

• If we are in the regime where p <C pc, this means that the quasiparticle has a speed much smaller than Cs, so 
that this is the Galilean limit of the relativistic regime. We are then dealing with a Galilean quasi-particle. The 
rest mass Mc^ provides the usual constant shift of the Gahlean energy £ ~ -Mc^ + ^p'^/M. 

II. GRAVITATIONAL DYNAMICS 

Since we have described the physics of the system in the case of homogeneous condensate ip, we can now pass 
to study the inhomogeneous one, and hence the emergence of a gravitational dynamics. We are going to consider 
condensates which are nearly, but not exactly, homogeneous: this will correspond to the case of weak gravitational 
field. 

In a curved spacetime, the identification of the Newtonian gravitational potential requires a non-relativistic limit 
of the geodesic equation in a weak gravitational field [14|. For instance, in asymptotically fiat spacetimes, there 
is a coordinate system such that the metric, in the asymptotic region, can be written as g^i, = "q^n + /i^^, where 
hf^i, encodes the deviation from fiatness, i.e. the gravitational field. In this regime, it is possible to show that the 
Newtonian gravitational field is identified with the component /iqo : 

$^(x) = -i/ioo(x) (31) 

It is well known that in the context of standard BEG (i.e. dealing with the non-linear equation ([3])), the quasi-particles 
travel in an emergent metric ds^ determined in terms of the homogenous condensate ij) ■ 



ds2 ^ J!^ [_ (^2 _ ^2-) ^^2 _ 2v,dtdyi' + JydxMxJ] , (32) 



where m is the mass of the atoms and Cg and v depend on the properties of the condensate ip = ^/rQe^^ , through 

unc ^ 1 - 

Cc — , V — — VO. 

m 771 

Considering that the condensate is homogenous, the density and velocity profiles become constant, i.e. respectively 
"^c — noo, V = ^oo- With the coordinate transformation, 

dT = dt, dX' = dx^ - v'^dt, (33) 

the line element ([32]) is rewritten as: 

dsl = -cldT^ + dX^. (34) 

The condition of asymptotic flatness for spacetimes can be translated with the condition of asymptotic homogeneity 
for the condensate. We require then that only in a small region of space, in the bulk, the condensate deviates from 
perfect homogeneity. 

We consider therefore some small deviation from the asymptotic values of the velocity and of the density: 

JT-c = "■00(1 + 2u(x)), v = Voo+w{x), with M <C 1, W'tiv. (35) 

This implies in particular a rescaling of the speed of sound. 

cL(l + 27.(x)). 
rn 

The acoustic line element (|32]) becomes then 



2 hiTlQ 2 



ds^ = J!^!!^ (_(c2 _ v^)dt^ - 2v,diedt + S,,diedxi) , (36) 

where we have introduced a constant prefactor mcoo/nao in order to have the conformal factor asymptotically nor- 
malized to one. Using ([55]) . together with the coordinate change (|33p . the acoustic line element has the form: 

ds^ = ds^ - iu{X)cl^dT'^ ~ 2w,{X)dTdX' + u{X)d^jdX'dX\ (37) 

at first order in u, Wi. Consequently, we see that 

h^oiX) = -3cLu(X), (38) 

so that the gravitational field is encoded in the density perturbation of the condensate wave-function ip, 

•^NiX) = \cl,u{X), (39) 

while it is independent from velocity perturbations. 

In light of this discussion, we shall discuss situations in which the condensate wavefunction has a constant phase, 
while its modulus slightly deviates from perfect homogeneity. In order to stay as close as possible to the case of the 
standard BEC analysis we have just presented, it is convenient to introduce the parametrization: 

V'-(^) (l + «(x)), (40) 

where u(x) is dimensionless and it is assumed to be very small. In practice, we will assume that it is associated with 
a localized inhomogeneity of the condensate. At infinity (i.e. at the boundary) we ask that u — > 0. 

A. Identifying the gravitational potential for the quasi-particles 



First of all we want to see if there is a term in the equation of motion for quasi-particles (J16p . which can be identified 
as an external potential term. This will allow us to identify the gravitational degrees of freedom. We have seen in the 
previous section that they should be encoded in the condensate wave- function -0. We note however that the dynamics 
of if) is essentially non relativistic (c.f. (|12p ). We are therefore looking for a Newtonian type of gravity. To identify the 



Newton potential, we therefore need to look at the Galilean regime for matter, i.e. in the regime where the momentum 
p of the quasi-particle is such that p <^ pc (c-f- section [TCt . 

Note that the usual BEG construction given in ^ gives rise to massless particles. However these latter cannot be 
handled in the context of Galilean mechanics. It is then not clear at all how one can identify a Poisson like equation 
for gravity in this case. To solve this issue, we have introduced Hx in order to generate a non zero mass for the 
quasi-particle. We can therefore now consider the quasi-particle as a possible source to the Newtonian gravitational 
potential. 

To do identify the Newton potential, we repeat the diagonalization of the Hamiltonian in (|16p for the field %, now 
including the fluctuations of the condensate. In this case, the diagonalization procedure is more involved: we have to 
deal with the non-commuting operators V^ and u. We can not perform it in an exact way. However, we are interested 
in the Galilean regime for the quasi-particle spectrum, when pc ^ P- It is then a reasonable approximation to neglect 
all the terms involving the commutators [p^/2to,m(x)], which are largely suppressed (with respect to the other terms 
appearing in the equations) by the mass of the atoms and from the smallness of m(x). 

With these simplifying assumptions, the Hamiltonian for the quasi-particles in the non-relativistic limit is 

H . «A^c^_^!Z! + 2l^i±M^i±i^^,(x) (41) 

where the mass of the quasi-particle M. and for the speed of sound Cs are given in (|26p . We first recognize the constant 
shift M.c^s of the energy due to the rest mass in the Galilean regime. This term is not affecting the discussion in any 
ways and can be eliminated. The term proportional to m(x) can be clearly interpreted as an external potential. If we 
want to identify it with the gravitational potential $grav, we need to have 

2^^^ J^2 ^"W = -^^grav ^ $grav(x) = ^ ^ ^m(x), (42) 

where Ai is the mass of the quasi-particles. Note that this identification is formal, and relies on the way in which the 
gravitational potential enters the Schroedinger equation for a non-relativistic quantum particle. We should always 
work with u: our definition of 'I'grav is dictated from the analogy we want to make with Newtonian gravity. For 
instance, we see that this definition becomes singular when we deal with massless quasi-particles, i.e. when A — > 0. 

B. (Modified) Poisson equation 

Now that we have identified a candidate for the Newton potential ^giav from the quasi-particles dynamics, we need 
to check that it satisfies some sort of Poisson equation. Since the gravitational potential is deduced from V' - as small 
deviations from perfect homogeneity (c.f. (|40|) ) - the Poisson equation should be deduced from the BdG equation (fT2|) . 
With the natural assumption that the potential is reacting instantaneously to the change of distribution of matter, 
we can neglect the time derivative and (|12p becomes 

^V' - 2(^1 + A)^ uix) = 2« (^n(x) + im(x)^ (43) 

We have seen in section II Al that the terms m(x) and n(x) are functions of the atoms x outside the condensate and 
therefore of the quasi-particle (p, through the Bogoloubov transformation (I20p . They can be therefore interpreted as 
the source in the (modified) Poisson equation. We examine now different types of source: either localized particles or 
plane- waves. 

1. Localized sources 

The most natural source to consider for the Poisson equation is a single quasi-particle </> at a given position xq. 
However, point-like particles give rise to divergences. We consider therefore a quasi-particle which is localized around 
the point xq, with a non-zero spread to regularize these divergencies. We consider a quasi-particle in a state of the 
form: 

ICxo)- /d3xCxo(x)<^Ux)|r!), with /d^x|Cxo(x)P = l^(CxolCxo> = l- (44) 



^xo encodes the spreading of the particle around xq since 

(C.o|(^nx)<^(x)|Cxo> = ICxo(x)p. (45) 

We can now determine the value for the anomalous mass m and anomalous densities n when the quasi-particle is in 
the state |Cxo)- An explicit calculation, given in the appendix [b1 gives 



n(x) 



2 2 

y'dV(x-z)Cxo(z) + y'd3zg(x-z)Cxo(z) +:^E/3'(k), (46) 

k 

m(x) - 2 (^1 d=^zi5(x - zi)C (zi)) (^J d'z.fix - z2)Cxo (^2)) + ^J2 "(k)/3(k), (47) 

where we have introduced the functions /, g depending on the Bogoliubov coefficients a and P 

/(x) = 1 E a(k)e^^-, 5(x) = -^ ^ /3(k)e-*-. (48) 

k k 

The quantities nn and mo with 

"" = ^ E /3'(k), nx,, = 1 E "(k)/3(k), (49) 

k k 

are vacuum contributions independent from the presence of actual quasi-particles. They are related to the inequiva- 
lence of the particle and quasi-particle vacua, and it can be easily seen that: 

nn = {n\xHx)x{x)\n), mo = (f^|x(x)x(x)|r!). (50) 

The functions /, g encode the fact that quasi-particles are collective degrees of freedom and therefore intrinsically 
some non-local objects. This non-locality is precisely encoded in the Bogoliubov transformation (PO)) . Quasi-particles 
and atoms (i.e. local particles) coincide only if we have a(k) — l,/3(k) = 0, and therefore /(x) = <5'^(x), while 
g{x) = 0. Since this is not the case, the anomalous mass and the anomalous density will show an intrinsic non- 
locality. The spreading characterized by |Cxo) encodes some extra non-local effect, i.e. the quasi-particle is in some 
sense an extended object. It was however introduced by hand, for a regularization purpose and therefore is not 
fundamental as the non-locality introduced by the Bogoliubov transformation. 
The equation (|43)) becomes then: 

— V^ - 2(/i + A) j u{x) = 2k (n{x) + ^m(x) j + 2k (nn + ^mnj , (51) 

where we have introduced the quantities 

n(x) =: n(x)-no, m(x) = m(x) - mn, (52) 

which represent the contribution of actual quasi-particles to the anomalous density and anomalous mass, respectively. 
By dimensional analysis, the terms n, m have the dimensions of a number density. Since in Newtonian gravity the 
source for the gravitational field is a mass density, we introduce the mass density distribution: 

Pmattcr(x) = M ( n(x) + ^"^(x) ) ■ (53) 



With this definition, we can rewrite ([5T|) as an equation for the field <&, 



grav ■ 



V^ - -^ ) $g,av = 47rG'^>„,attcr + A, (54) 



where we have defined 



Woc_ At(/i + 4A)(M + 2A)^ 2k(m_4A)(a^+2A) 1 

^^ = 4.n2™A3/2(^ + A)V2' ^ ft^A ("^^ + 2"^^)' (^^) 

L^^ , , ,, . (56) 

4m(^ + A) ^ ' 
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This particular choice of notation is motivated by the comparison of ((54)) with the Newtonian hmit of Einstein 
equations with a cosmological constant (see, for instance, Eq. (9) of reference [I^l). For this reason, we can identify 
these three quantities as the analogous of the Newton constant, the analogous of the cosmological constant and a 
length scale which represents the range of the interaction, as we are going to discuss below. 

To get a better grasp of the physics of the modified Poisson equation ([54)) , we can look at its solution for a given 
distribution of quasi-particle pmattcr- 

As it is well known, a solution for the equation 

V^ " Jl) *W = 4^G^'A^'5'(x - z), (57) 



is given by the Yukawa potential 

c&Y (x; z) = ^^^^ii^, . (58) 

On the other hand, a solution for the equation 



G^^^Me-^""-^^/^ 



V^ - j^j <i>(x) = A, (59) 

is just given by the constant solution 

$A = -i^A. (60) 

The linearity of equation ()M)) allows us to use these results to write down a solution for a generic distribution of 
matter (i.e. quasi-particles) as 

$grav(x) ^ / /9mattcr(z)$Y(x; z)d^Z + $A. (61) 

Solutions of (|54)) are therefore constructed from the Yukawa potential smeared out due to the non-locality of the quasi- 
particle (with an extra global shift due to the cosmological constant). The Yukawa potential is typically encoding 
some short range interaction, characterized by the scale L which is simply related to the healing length (fTS)) . 



i'-y- (62) 

Although this a very short range for the gravitational interaction, this outcome should not come as a surprise. In 
fact, the healing length (c.f. (fT5| ) characterizes the typical length over which a condensate can adjust to density 
gradients. Since density inhomogeneities encode the gravitational interaction, one should expect them to be damped 
over a distance of the order the healing length. 

In the context of relativistic field theory, the short interaction scale for gravity would be translated in a massive 
graviton, with mass given by 

, r9 ^^ . M + -^ 2 
M = = 4— m 

We can then compare the masses of the quasi-particles M, graviton Mgrav and atoms m, 

< 7W < m < V2m < Mgrav < 2m, (63) 

which shows the hierarchy of the energy scales present in this system. We notice that the graviton is then always 
more massive than the quasi-particles, and that this interaction is of very short range, since the ^ is much shorter 
than the acoustic Compton length^ of the quasi-particles. In particular, we cannot tune the parameters of the system 
in such a way to make Mgrav arbitrarily small, in order to be closer to reality. 



We are using Cgo^nd instead of cught to define all these scales. We have to use the natural units for an hypothetical phononic observer. 
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2. Plane waves 

While a quasi-particlc localized in a given point in space is certainly the most natural source for gravity from the 
Newtonian perspective, it is interesting also to see what happens when instead we consider quasi-particles with a 
definite momentum p = ?ik. Let us focus first on the special case of a 1-particle state with momentum p, that is 
|p) = b\\^)- The anomalous mass and the anomalous densities become then 

n(x) = (f}|Stx^(x)x(x)6k|fi), m(x) = {^\hlx{^)x{^)k\n) . (64) 

To express them in terms of quasi-particles, we need to perform the Bogoliubov transformations (|^D|) . As we recalled 
in the previous section, we can not specify these transformations exactly, due to the presence of the potential u. 
However, the corrections to the Bogoliubov coefficients a{k), (3{k), evaluated in the case u = in ((23|) . provide 
corrections to the expressions of n and tn which are relevant only beyond the linear order in u, which means beyond 
the approximation we are using. We can then safely neglect these corrections. Using the Bogoliubov transformation, 
we obtain explicitely 

, , a2(fc)+/?2(fc) ,a(fc)/3(fc) 

nk{x)= ^ ' ^"^ ^ ^ +no, mfc(x) = 2 ^ '^ ' +mn. 

where we recognize the contribution of the vacuum no , mn . 

The generalization to the case of states containing a definite number of quasi-particles with a given momentum 
follows in the same way. For these states denoted as \n{ki), ...,7i(fc„)), one obtains: 

n,(x) ^^n(fe.) "'^'-^;^'^'-^ n., m.(.) = 2 5:^(^0^^^^^^ + mo. (65) 

i i 

In all these expressions, besides the Bogoliubov coefficients, we recognize the terms n(ki)/V, which are the number 
densities of quasi-particles in a given eigenstate of momentum. These number densities, however, are not giving 
immediately the source term to be inserted into (|M|) . since they are weighted by the Bogoliubov coefficients. This is 
the representation, in momentum space, of the non-locality we have discussed in position space. 
We rewrite (j54p with a source term made by a single particle with a given momentum p = fik. 

V^ - — j $grav =- 47rG'Ar(fc)p„,attcr + A, (66) 

where we have pmattcr = M/V since we have just one particle of mass A4, while A, L are defined as in ([55l) - (f56|) . We 
encode the effect of the Bogoliubov coefficients, and hence of non-locality, in the "running" Newton constant 

Gjv(fc) = (a2(fc) + /32(fc) -I- a(fc)/3(fc)) G%^. (67) 

The discussion of the solution to this equation is even simpler than the localized state, given that the source term is 
just constant. Consequently, <l>grav = const is a solution as in the case with purely vacuum contribution. For what 
concerns the physical effects of this kind of gravitational field, we have to plug this constant solution into (HT|) : this 
amounts just to a shift of the energy, which is, however, momentum dependent, leading to observable relative energy 
shifts if different momenta are considered. 

Conclusions and Remarks 

In an analogue gravity model based on a BEC system, the degrees of freedom are separated into the atoms that 
condense and the one which do not. Quasi-particles are then collective degrees of freedom constructed from the 
un-condensed atoms. The dynamics of the quasi-particles is encoded, in a given regime, as particles propagating in a 
curved spacetime metric, which is characterized by the the density Uc and the velocity profile v of the condensate. In 
this sense, we can expect that gravitational degrees of freedom are encoded in the condensate. Dynamics of the latter 
is encoded in the BdG equation p^ . which is essentially Galilean. Hence, we can not expect to recover the Einstein 
equations in this context [lO[. Nonetheless, one can still try to interpret (|12p as some sort of Poisson equation for 
some type of Newtonian gravity. 

However, quasi-particles are massless in usual BEC systems and hence they cannot be considered as sources for 
the gravitational field in the Poisson equation. We introduced therefore a new term Hx in the dynamics of the BEC 
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which softly breaks the U{1) symmetry and consequently, as we showed in section II Bl generates a mass gap for 
the quasi-particles. We showed explicitly that the presence of this small symmetry breaking term does not prevent a 
condensation from happening and still allows a mean field description (which is sufficiently accurate for our purposes). 
Then, following the usual General Relativity argument, we have argued, in section pT)) . that the Newtonian potential 
$7V has to be related to small inhomogeneities in the condensate density (while perturbations in the velocity profile do 
not contribute at first order as gravitational degrees of freedom). This conjecture, based on the analysis of a standard 
BEG system, was then confirmed by a specific analysis of the modified BEG dynamics for an almost homogenous 
condensate. 

The end point of this investigation can be then summarized in the following two equations 

F = Md= -TWV^grav, (68) 

V2 - -l") $g,.av - AttGnP + A, (69) 

where M is the mass of the quasi-particle acquired via the soft C/(l) symmetry breaking induced by (O, L is 
proportional to the healing length, A plays the role of the cosmological constant and Gn is an effective coupling 
constant that depends on the condensate microphysics and the form of the matter source. 

For what regards the latter we have considered two cases: a localized quasi-particle state and a set of plane waves. 
In the first case the analogue Newton constant is indeed momentum and position independent and the solution of the 
modified Poisson equation ([M|) . has the form of a smeared Yukawa potential. The smearing is due to the fact that 
quasi-particles are intrinsically non-local objects, being collective degrees of freedom. When considering plane-waves 
as sources, we have instead that, due to the momentum dependence of the Bogoliubov transformation, Gn is running 
with the momentum and the solution for the gravitational potential is a constant (albeit a different one for different 
momenta). One should however be careful: while it is common in quantum field theory (QFT) to encounter the 
notion of running coupling constants, the origin of the running here is rather peculiar. Indeed, in QFT the running 
is due to quantum corrections to the tree level/classical action, here the running is due to the inequivalence between 
the ground state of the Fock spaces of atoms and quasi-particles. Paraphrasing what has been done in the context 
of emergent geometry, where the notion of "rainbow geometry" has been introduced, we could speak about "rainbow 
dynamics" . 

We have also obtained naturally a cosmological constant in the model: vacuum gravitates. It is induced by the 
terms (ri|x^xl^)j (^IxxI^)j where 51 is the state with no quasi-particles. It is entirely due to the (unavoidable) 
inequivalence between the quasi-particle vacuum and the particle vacuum and cannot be put to zero just tuning the 
parameters. It represents an interesting alternative to known mechanisms to generate a cosmological constant (see 
also [16| for similar ideas about the nature of the vacuum energy in condensed matter systems) . 

In conclusion, BEG as an analogue model for gravity presents therefore many differences with a realistic gravity 
theory. One should not be deceived by this result, as it would have been a preposterous expectation to recover 
Einstein General Relativity in a Bose Einstein Gondensate. The model is however interesting per se as it still encodes 
a modified Poisson equation and hence provides new insights on the possible origin of the cosmological and Newton 
constants in emergent gravity scenarios. 

As a future development of this investigation, it might be interesting to analyze a 2-BEG model [Ij]: in fact in 
this case one could treat a multi-particle system whose richness could allow a closer mimicking of Newtonian gravity. 
However, the fact that emergent gravity has to be Newtonian in a BEG based analogue model seems to be unavoidable 
since the gravitational potential depends on the condensate, which is typically described by non-relativistic equations. 
A possible way to avoid this issue is either to consider relativistic BEG [l^ (however in this case we would still 
expect to get only some type of scalar gravity), or to change completely paradigm and identify gravity not in the 
condensate but among the perturbations around the condensate (see for example [IS])- We leave these ideas for 
further investigations. 
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APPENDIX A: THE FLUID DESCRIPTION 

The Gross-Pitaevskii (GP) equation deseribing a BEG admits an interesting fluid interpretation, through the 
Madelung representation. We are considering the GP equation given in (|13p 

ot 2m 

and we want to use the Madelung representation for the complex field ip: 

^ = y^,e-''/\ (A2) 

When replacing this into the GP equation, dividing by the phase and splitting the resulting expression into the real 
and imaginary parts we obtain two equations: 

Tic + V • [ncV) = -T"-c ^i'^ ( "T ) ' y^^) 

Q = Vijuantum + -^^'^'^ - /i - A COS ( — I - KTlc, (A4) 



2 

where we have introduced the velocity field v — —V9/h, and 



2 



1 h ., 

V quantum ^ ^= 7^ ^ v'^o l-^&j 

^/r^r 2?77, 



is the familiar quantum potential term. These two equations, in the case A = 0, have a nice interpretation as the 
continuity equation and the Euler equation for a perfect fluid. The continuity equation, in particular, is just the 
statement about the conservation of the Noether current associated with the f/(l) invariance of the system. On 
the other hand, when A ^ the C/(l) invariance is broken, and the number operator is no more conserved by the 
Hamiltonian evolution. 

It is interesting to see what happens when we consider the case of homogeneous condensates, dyjic = d^^v^ = 0. 
From the first equation we get: 

sin f — j =04^61 = y/i, ; e Z. (A6) 

This result implies that not only v is constant, but that actually vanishes. Inserting this result in the second equation 
we obtain: 

^^^^ + cos(/.)A^ (A7) 



From the analysis of the quasi-particle dynamics in section II Bl we have seen that A < corresponds to a negative 
mass square, i.e. tachyonic behavior. Since cos(/7r) = — 1 would be equivalent to changing the sign of A, without 
repeating the analysis of section HI we see that cos(Z7r) = 1 is required for stability (no tachyons). 

APPENDIX B: SOURCE TERM 

In this section we provide the details of the calculation of the source term for the Poisson equation corresponding 
to a localized source (c.f. IIIBT]) . We have to evaluate the expressions: 

n(x) = (Cxolx^(x)x(x)|ao), m(x) = (ajx(x)x(x)|ao), (Bl) 

where 

1^0)= /d=^zC.o(x)0^(z)|f^)- (B2) 
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Let us describe it for n, since m can be evaluated following the same steps. First, one has to write explicitly n in 
terms of the field operators: 

d3zid3z2C(zi)Cxo(z2)(r!|0(zi)x^(x)x(x)^t(z2)|r!). (B3) 

Let us evaluate then the expectation value inside the integral. To do this, it is necessary to replace the expansion of 
the field operators in plane waves, and then to use the Bogoliubov transformation: 



(f^l0(zi)xt(x)x(x)^t(z2)|f7) = -l J2 ^*''^ 



-ik-x ik'-x -ih'-Z2 „ 

2 . - e e e x 

k.k'.h.h' 



x(f7|6h(a(k)fej, + /3(k)Lk)(a(k')V +/3(k')foLk')^Ll^)- (B4) 

It is easy to see that, in this last expression, there are only two non- vanishing terms 

{v) = (f]|foh(a(k)Sj^ +/3(k)Lk)(a(k')6k' +/3(k')foLk')^Ll^> = «(k)a(k')(r!|6h&I&k'&i|0) + 

+/3(k)/3(k')(f^|&h^-k6Lki|j^>- (B5) 

Using the algebra of the operators 6, w , it is easy to see that the expression reduces to: 

(v) = a(k)a(k')<5h,kV.k' + /3(k)/3(k') ((5k,k"5h.h' + <5h,-k'<5k,-h') • (B6) 

Consequently, 

n(x) = A{x) + B(x) + C(x) , (B7) 



where 



^W = T72 /d'zid'z2 J2 e*-^^e-*"^e*'"^e-*'-^^Co(zi)C.o(z2)a(k)a(k')<5h,k<5h',k', (B8) 

•^ h hi h hi 



k,k',h,h 



Bix) ^ ^ f d'z,d'z2 J2 e*-^^e-*-V^'-^e-*'-^X*„(zi)C.o(z2)/3(k)/3(k')4,k'4,h' 

■^ h Irl V, hi 



(B9) 



k.k'.h.h 



^W = T72 /d'^id^zs J2 e*-^^e-*-^e*'"^e-*'-^Co(zi)C.o(z2)/3(k)/3(k')<Jh,-k'4,-h'. (BIO) 

•' k,k',h,h' 

To manipulate these expression, it is useful to recall the representation of the Dirac delta in a Fourier series: 

'5'(xi-x2) = -i5]e-*(-^--^), (Bll) 

k 

and that the distribution (^^o is normalized, 

y'd3x|C.o(x)p^l. (B12) 

After straightforward passages we obtain 



A(x) = 



2 

d3z/(x-z)C.„(z) 



(B13) 



^W = ^E/^'(k), (B14) 

k 
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and, finally: 



n(x) 



C(x) = 



d^zf{x-z)C^„{z) 



d^z5(x-z)C^o(z) 



d^z5(x-z)Ci;o(z) 



where we have introduced the functions: 

/(x) = 1 ^ aik)e^^-\ 5(x) = ^J2 ^(k)^ 



^E/5^(k) 



V 



— ik-yi 



(B15) 



(B16) 



(B17) 



Notice that, as a consequence of a(k) = q;(— k), /3(k) = /?(— k) and of the fact that these coefficients can be chosen to 
be real, the functions f,g are real functions. 

Applying the same procedure to the term m(x) we obtain 

m(x) = 2(^|d3zig(x-zi)C„(zi)) (^|d3zi/(x - Z2)C.„(z2)) + -^ ^a(k)/3(k). (B18) 
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